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Abstract 

We consider general second order uniformly elliptic operators subject 
to homogeneous boundary conditions on open sets <^(r2) parametrized by 
Lipschitz homeomorphisms (p defined on a fixed reference domain 0. Given 
two open sets we estimate the variation of resolvents, eigenvalues 

and eigenfunctions via the Sobolev norm ll'?^* — for finite values of p, 
under natural summability conditions on eigenfunctions and their gradients. 
We prove that such conditions are satisfied for a wide class of operators 
and open sets, including open sets with Lipschitz continuous boundaries. 
We apply these estimates to control the variation of the eigenvalues and 
eigenfunctions via the measure of the symmetric difference of the open sets. 
We also discuss an application to the stability of solutions to the Poisson 
problem. 
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1 Introduction 

This paper is devoted to the proof of stabihty estimates for the non-negative self- 
adjoint operator 

i,j=l * ^ 

subject to homogeneous boundary conditions, upon variation of the open set fl in 
R^. Here A, are fixed bounded measurable real-valued functions defined in R 
satisfying = Aji and a uniform ellipticity condition. 
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The focus is on explicit quantitative estimates for the variation of the resolvents, 
eigenvalues and eigenfunctions of L on a class of open sets diffeomorphic to Q. 

In the first part of the paper we consider two diffeomorphisms 0, of onto 
4>{fl), respectively, and we compare the resolvents, eigenvalues and eigen- 

functions of L on the open set with those of L on 0(f2). The main aim is to 
provide stability estimates via ||0 — 0||vi/i.p(n) for finite values of p. These estimates 
are applied in the last part of the paper where we take (p = Id and, given a de- 
formation Cl of Q, we construct a special diffeomorphism (p representing Cl in the 
form Q = (f){Q) and obtain stability estimates in terms of the Lebesgue measure 
|i7 A f2| of the symmetric difference of Q and Q. 

Our method allows us to treat the general case of the mixed homogeneous 
Dirichlet-Neumann boundary conditions 

M = on r, and ^ AjjT^^i = on 91] \ T, (1.2) 

where F C dQ. To our knowledge, our results are new also for Dirichlet, and for 
Neumann boundary conditions. 

There is a vast literature concerning domain perturbation problems, see for in- 
stance the extensive monograph by Henry [15] and the survey paper by Hale [14] . 



In particular, the problem of finding explicit quantitative estimates for the varia- 
tion of the eigenvalues of elliptic operators has been considered in Burenkov and 
Davies [3], Burenkov and Lamberti [U El [6], Burenkov and Lanza de Cristoro- 
foris [S], Davies [SI [ID], Lamberti and Lanza de Cristoforis [ISlilZl, and Pang [2U] : 
see Burenkov, Lamberti and Lanza de Cristoforis [7] for a survey on the results 
of these papers. However, less attention has been devoted to the problem of find- 
ing explicit estimates for the variation of the eigenfunctions. With regard to this, 
we mention the estimate in [20] concerning the first eigenfunction of the Dirichlet 
Laplacian and the estimates in [THlilT] concerning the variation of the eigenprojec- 
tors of the Dirichlet and Neumann Laplacian. In particular, in [TBlfTT] the variation 
of the eigenvalues and eigenprojectors of the Laplace operator was estimated via 
II V0 — V0||L°°(n) under minimal assumptions on the regularity of fl, (p and (p. 

In this paper, we consider the same class of transformations (p, (p as in [TBI [12] 
(0, bi-Lipschitz homeomorphisms) and by making stronger regularity assump- 
tions on 0(fi), 0(fi) we estimate the variation of the resolvents, eigenvalues, eigen- 
projectors and eigenfunctions of L via the measure of vicinity 

(5p(0, 0) := II V0 - V0||LP(n) + II A o - A o (P\\Lp{n) (1.3) 

for any p e]po, oc], where A = (Aj)ij=i,...,Ar is the matrix of the coefficients. 

Here po ^ 2 is a constant depending on the regularity assumptions. The best 
Po that we obtain is po = N which corresponds to the highest degree of regularity 
(see Remark (TUI), whilst the case po = oo corresponds to the lowest degree of 
regularity in which case only the exponent p = oo can be considered. The regularity 
assumptions are expressed in terms of summability properties of the eigenfunctions 
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and their gradients, see Definition HI Observe that if the coefficients Aij of the 
operator L are Lipschitz continuous, then 5p(0, 0) does not exceed a constant 
independent of 0,0 multiphed by the Sobolev norm ||0 — (p\\w^'P{n)- Moreover if 
the coefficients A^j are constant then the second summand in the right-hand side 
of fll.3p vanishes. 

More precisely, we prove stability estimates for the resolvents in the Schatten 
classes (Theorem [8]), stability estimates for eigenvalues (Theorem [T3|) . eigenpro- 
jectors (Theorem [T6|) . and eigenfunctions (Theorem [201) . In the Appendix we also 
consider an application to the Poisson problem (we refer to Savare and Schim- 
perna [22] for stability estimates for the solutions to the Poisson problem in the 
case of Dirichlet boundary conditions obtained by a different approach). To prove 
the resolvent stability estimates in the Schatten classes we follow the method de- 
veloped in Barbatis [H [2] . 

In Section [7] we apply our general results and, for a given deformation Q of Q, 
we prove stability estimates in terms of \Q A Cl\. This is done in two cases: the 
case in which Q is obtained by a localized deformation of the boundary of Q and 
the case in which is a deformation of Q along its normals. We also require that 
the deformation F of F is induced by the deformation of Q (see conditions (17.31) 
and (I7.14P ). In these cases, similarly to [5], we can construct special bi-Lipschitz 
transformations (p : Q ^ Q such that 0(F) = F and 

||V0-/||lp(q) <c|fi Af]|^/^ (1.4) 

Observe that using finite values of p is essential, since in the case p = oo the 
exponent in the right-hand side of (11.41) vanishes. 

Let us describe these results in the regular case in which Q, Q are of class 
C^'^ and F, F are connected components of the corresponding boundaries. In 
Theorems [29], [32] we prove that for any r > N there exists a constant Ci > such 
that 

1 1 



E 

vn=l 



A„ + 1 A„, + 1 



<ci\nAn\-, (1.5) 



if |f2 A r2| < c|f^. Here A„, A„ are the eigenvalues of the operators (11. ip corre- 
sponding to the domains Q, Q and the associated portions of the boundaries F, F 
respectively. Moreover, for a fixed Q and for any r > N there exists C2 > such 
that if A„ = . . . = Xn+m-i is an eigenvalue of multiplicity m then for any choice 
of orthonormal eigenfunctions ipn, • • • , V'n+m-i corresponding to A„, . . . , Xn+m-i, 
there exist orthonormal eigenfunctions ipn, ■ ■ ■ , ipn+m-i corresponding to A„, . . . , 
Xn+m-i such thal0 

ll^fc - V^fc|lL2(nun) < C2|l^ A (1.6) 

for all A; = n, . . . , n + m — 1, provided that |fi A f2| < ^. Here it is understood 
that the eigenfunctions are extended by zero outside their domains of definition. 



^Note that, for a fixed and variable (l, one first chooses eigenfunctions in fl and then finds 
eigenfunctions in il, while the opposite is clearly not possible. 
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In the general case of open sets fi, with Lipschitz continuous boundaries and 
r, r with Lipschitz continuous boundaries in dQ, dQ our statements still hold for 
a possibly worse range of exponents (see Theorems [291 E2D- 

We emphasize that, in the spirit of [121 [IZ], in this paper we never assume 
that the transformation belongs to a family of transformations (pt depending 
analytically on one scalar parameter t, as often done in the literature (see e.g., [T3] 
and [H] for references). In that case, one can use proper methods of bifurcation 
theory in order to prove existence of branches of eigenvalues and eigenfunctions 
depending analytically on t. In this paper is an arbitrary perturbation of and 
this requires a totally different approach. 

The paper is organized as follows: in Section [2] we describe the general setting; 
in Section [3] we describe our perturbation problem; in Section H] we prove stability 
estimates for the resolvents and the eigenvalues; in Section |5] we prove stability 
estimates for the eigenprojectors and eigenfunctions; in Section [6] we give sufficient 
conditions providing the required regularity of the eigenfunctions; in Section [7] we 
prove stability estimates via the Lebesgue measure of the symmetric difference of 
sets; in the Appendix we briefly discuss the Poisson problem. 



2 General setting 

Let n be a domain in of flnite measure. We consider a family of open sets 
[Q) in parametrized by Lipschitz homeomorphisms of onto 0(f2). Namely, 
following [16], we consider the family of transformations 

$(f2) := |0 G (L^'°°(n))^ : the continuous representative of 

is injective, essinf | det V0| > o| , (2.1) 

where L^'°^(i7) denotes the space of the functions in Lj^^{Q) which have weak 
derivatives of first order in L°° (Q). Observe that if G $(f2) then is Lipschitz 
continuous with respect to the geodesic distance in Q. 

Note that if G $(i^) then 0(f2) is open, is a homeomorphism of Q onto 
0(f2) and the inverse vector-function (p^'^^ of belongs to $ (0(f2)). Moreover, any 
transformation G $(^^) allows changing variables in integrals. Accordingly, the 
operator from L^(0(i7)) to L^(fi) defined by 

is a linear homeomorphism which restricts to a linear homeomorphism of the space 
Vri'2(0(fi)) onto Vri'2(fi), and of V^o'^(0(^])) onto W^^'^ip). Furthermore, \/{vo 
0) = Vt;(0) V0 for all v G VT^'^ (0 (f2)). Observe that if G <^>(^]) then the measure 
of 0(fi) is finite. See [IH] for details. 
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Let A = (Aij)jj=i,...,Ar be a real symmetric matrix-valued measurable function 
defined on such that for some 9 > 

N 

e-'\^\'<Y,A,,{x)U,<0\^\\ (2.2) 

for all X, ^ G R^. Note that (Q implies that Aij e L°°(R^) for all i,j = 1,...,N. 
Let (p G and let W be a closed subspace of W^'^{(f){Q)) containing 

We consider the non-negative self-adjoint operator L on 
canonically associated with the sesquilinear form Ql given by 

Dom(Qi) = W, Ql{vi,V2)= V^A.^^dy, vi,V2eW. (2.3) 



Recall that v G Dom(L) if and only if f G W and there exists / G L^(0(n)) such 
that 

QL(t^,V^) = (/,V^)L2(^(f,)), (2.4) 

for all if) G W, in which case Lv = f (see e.g., Davies [H]). The space W is 
related to the boundary conditions. For example if W = WQ^'^(0(i7)) (respectively, 
W = Vr^'^(0(f2))) then the operator L satisfies homogeneous Dirichlet (respec- 
tively, homogeneous Neumann) boundary conditions. 

We also consider the operator H on L'^{Q) obtained by puUing-back L to L^(fi) 
as follows. Let v G W^''^{(f){Q)) be given and let u = v o (j). Observe that 



/ \v\'^dy = / det V0I (ix . 

J6(n) Jq 



'(j}{n) Jn 
Moreover a simple computation shows that 

f . 8v 8v ^ f ^ 8u 8u . ^ ^ , I 7 

where a = {aij)ij=i^,,,^N is the symmetric matrix-valued function defined on Q by 



Oil — ^ ^ ( ^ 



<92/r- dys 



r,s=l 

= ((V0)-iA(0)(V0)-*).,. (2.5) 

The operator H is defined as the non-negative self-adjoint operator on the Hilbert 
space L^(n, | det V0| dx) associated with the sesquilinear form Qh given by 

r ^ 8u 8u 

Dom(QH) = Q[>V], Qh{ui,U2) = / Oij^^l det V0| rfx, mi, M2 G C0[yV]. 
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Formally, 

Alternatively, the operator H can be defined as 

H = C^LC^(-i). 

In particular H and L are unitarily equivalent and the operator H has compact 
resolvent if and only if L has compact resolvent. (Observe that the embedding 
W C L^(0(n)) is compact if and only if the embedding C(^[W] C L'^{fl) is compact.) 

We set g{x) := | det V0(a;)|, x E Q, and we denote by (-, the inner product 
in L'^{fl,gdx) and also in {L'^{fl, g dx))^ . 

Let T : L'^{Q, g dx) —>■ {L'^{Q, g dx))^ be the operator defined by 

Dom(T) = C^W], Tu = a^'^Vu, u G C^[W]. 
It is then easily seen that 

H = T*T, 

where the adjoint T* of T is understood with respect to the inner products of 
L\n,gdx) and {L\n,gdx))^. 



3 Perturbation of cj) 

In this paper we study the variation of the operator L defined by (12.31) upon 
variation of (p. Our estimates depend on essinf^ | det V0| and || V0||Loo(f^). Thus 
in order to obtain uniform estimates it is convenient to consider the families of 
transformations 

<l>r{ft) = 1^ e <^{Q) : r"^ < essinf |detV0| and ||V0||L-{n) < r| , 

for all r > 0, as in [16]. Here and below for a matrix-valued function B{x), x & Q, 
we set 1 1 -B 1 1 LP (n) = || where |-B(a;)| denotes the operator norm of the 

matrix B{x). 

Let0,0G <l>,(fi). Let >V and >V be closed subspaces of iyi'2(0((])), (0(1])) 
respectively, containing WQ'^(0(f2)), Wq"^ respectively. We use tildes to 
distinguish objects induced by 0, VV from those induced by (p, W. We consider the 
operators L and L defined on L^(0(f2)), L^(0(f2)) respectively, as in Section [2l 

In order to compare L and L we shall make a 'compatibility' assumption on 
the respective boundary conditions; namely, we shall assume that 

C4W]=CAW]. (3.1) 
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This means that Dom((5_f/) = Dom(Qj^), a property which is important in what 
follows. Clearly (13.11) holds if either L and L both satisfy homogeneous Dirich- 
let boundary conditions or they both satisfy homogeneous Neumann boundary 
conditions. 

We shall always assume that the spaces W, VV are compactly embedded in 
L^((/)(n)), L^(0(f2)) respectively, or equivalently that the space V := C<^[W] = 
C^[VV] is compactly embedded in L^(fi). 

Moreover, we require that the non-zero eigenvalues A„ of the Laplace operator 
associated in with the quadratic form \ Vu\'^dx, m G V, defined on V, 

satisfy the condition 

c*:=5^Ar<oo, (3.2) 

for some fixed a > 0. (This is in fact a very weak condition on the regularity of 
the set VL and the associated boundary conditions.) 

For brevity, we shall refer to assumption (A) as the following set of conditions 
which summarize the setting described above: 



(A) : \ V := C^[W] = C^[W] is compactly embedded in L\Q), 
y condition (13. 2p holds. 

Remark 1 We note that ifQ is a domain of class C^'^, i.e., Q is locally a subgraph 
of Lipschitz continuous functions, then inequality liS.S^) holds for any a > N/2 (see 
e.g., 13], Netrusov and Safarov ^19] and also Remark\2^ below). We also note that 
by the Min-Max Principle |77l p. 5] and by comparing with the Dirichlet Laplacian 
on a ball contained in VL, condition ^3. 2) does not hold for a < N/2 (no matter 
whether Q is regular or not). 

In order to compare L and L, we shall compare the respective pull-backs H 
and H. Since these act on different Hilbert spaces - L'^{Q, g dx) and L'^{Q,g dx) - 
we shall use the canonical unitary operator, 

w : L?{VL, g dx) — ^iF'iVL.gdx) , u wu , 

defined as the multiplication by the function w := g^^'^g~^^'^- We also introduce 
the multiplication operator S on (L^(il))^ by the symmetric matrix 

w-^a-'/^da-'/^ , (3.3) 

where the matrix a is defined by (12. 5 p and the matrix d is defined in the same way 
with replacing 0. If there is no ambiguity we shall denote the matrix (13. 3p also 
by 5. 

As it will be clear in the sequel, in order to compare H and H we shall need 
an auxiliary operator. Namely we shall consider the operator T*ST, which is the 
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non-negative self-adjoint operator in L'^{Q, g dx) canonically associated with the 
sesquihnear form 

/ (aVwi ■ Vu2)gdx , mi, M2 G V . 
Jn 

It is easily seen that the operator T*ST is the pull-back to Q via of the operator 

L := ^ L (3.4) 

defined on L'^{(f){Q)). Thus in the sequel we shall deal with the operators L, L and 
L and the respective pull-backs H, H and T*ST. We shall repeatedly use the fact 
that these are pairwise unitarily equivalent. 

We denote by A„[i?], n G N, the eigenvalues of a non- negative self-adjoint 
operator E with compact resolvent, arranged in non- decreasing order and repeated 
according to multiplicity, and by ipnlE], n G N, a corresponding orthonormal 
sequence of eigenf unctions. 

Lemma 2 Let (A) be satisfied. Then the operators L, L, L, H , H and T*ST have 
compact resolvents and the corresponding non-zero eigenvalues satisfy the inequality 

J2 A„[i5;]-" < cc*, (3.5) 

for E = L,L, L, H, H, T*ST, where c depends only on N, r, 9. 

Proof. We prove inequality (13. 5p only for E = T*ST, the other cases being similar. 
Observe that the Rayleigh quotient corresponding to T*ST is given by 

{T*STu,u)^ _ {STu.Tu)^ _ J^jaVu ■ Vuygdx 
{u,u)g {u,u)g J^\u\'^gdx 

Then inequality (13. 5p easily follows by observing that 

fiVu ■ Vu > ri|(V0)"Vnp > ^"V"2|Vup, 

|detV0| < Ar!|V0|^ (3.6) 
and using the Min-Max Principle [TT| p. 5]. □ 



4 Stability estimates for the resolvent and the 
eigenvalues 

The following lemma is based on the well-known commutation formula (14. 3 p (see 
Deift 
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Lemma 3 Let (A) be satisfied. Then for all^eC\ {a{H) U a{H) U a{T*ST)), 
{w-'Hw - 0"' -{H- 0"' = Ai + A2 + A3 + B, (4.1) 

where 

Ai = {l-w){wT*STw-0~\ 
A2 = w{wT*STw-0~\^-w), 
A3 = -^{T*ST-0~\w-w'^){wT*STw-0'^w, 
B = T*S^/\S^/^TT*S^/^ -0'^S^/\S-^ - I){TT* -^)-^T . 

Proof. It suffices to prove (14.11) for ,^ 7^ 0, since the case ^ = ^ <^{H) Ua{H) U 
a{T*ST) is then obtained by letting 1^ — 0. 

Recall that T*T = H. Similarly T*T = H, where we have emphasized the 
dependence of the adjoint operation on the specific inner-product used. In this re- 
spect we note that the two adjoints of an operator E are related by the conjugation 
relation E* = w'^E*w~'^. This will allow us to use only * and not *. 

Observe that 

fj = {h^/^wyd'^^V = ti;2(sV2v)*^-2~i/2y _ w^T*ST . (4.2) 
Therefore, by simple computations, we obtain 

= w-\H - O-'w ~ {H - 

= w-\w^T*ST-^)-^w-{T*T-^)-^ 

= w-\w^T*ST - ^)-^w - {T*ST - ^w-^)-^w-^ 

+ {T*ST - ^w-^)-^w-^ - {T*ST - ^w-Y^ 

+ {T*ST - ew-2)-i - {T*ST - ^)-^ 

+ {T*ST - ^)-^ - {T*T - 
= + + A3 + {{T*ST - - {T*T - 0"') . 

To compute the last term we use the commutation formula 

-i{E*E-C)-^ + E*{EE* -C)-^E = I (4.3) 

which holds for any closed and densely defined operator E, see [I2]. We write (14.31) 
first for E = T, then for E = S^^'^T, and then we subtract the two relations. After 
some simple calculations we obtain {T* ST — C,)~^ — {T*T — ^)~^ = B, as required. □ 

We now introduce a regularity property which will be important for our esti- 
mates. Sufficient conditions for its validity will be given in Section O 

Definition 4 Let U be an open set in and let E be a non-negative self-adjoint 
operator on L'^(U) with compact resolvent and Dom(£') C W^'^{U). We say that E 
satisfies property (P) if there exist go > 2, 7 > 0, C > such that the eigenf unctions 
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of E satisfy the following conditions: 

Unm\L^oiU)<CX^[E]^ (PI) 

and 

\\VME]\\L'>oiu)<CXn[E]'^+"2 (P2) 
for all n E N such that \n[E] ^ 0. 

Remark 5 It is known that iffl, Aij and T are sufficiently smooth then for the op- 
erator L in U.l\) . subject to the boundary conditions U.2\) . property (P) is satisfied 
with go = oo and 7 = N/4; see Theorem 24 and the proof of Theorem [2R 



By interpolation it follows that if conditions (PI) and (P2) are satisfied then 

10(1-^) 90(9-2)7 90(9-2) 1 I 90(9-2)7 

\\ME]\\lHU) < C''(«>-2)A„[E] «(^0-2) , \\Vijn[E]\\L.(U) < C'^(«0-2)A„[E]5+ 9(90-2) , 

(4.4) 

for all q E [2, go]- 

In the sequel we require that property (P) is satisfied by the operators H, H 
and T*ST which, according to the following lemma, is equivalent to requiring that 
property (P) is satisfied by the operators L, L and L respectively. 

Lemma 6 Let (A) be satisfied. Then the operators H , H and T*ST respectively, 
satisfy property (P) for some go > 2 and 7 > z/ and only if the operators L, L 
and L respectively, satisfy property (P) for the same go and 7. 

Let E' be a non-negative self-adjoint operator on a Hilbert space the spectrum 
of which consists of isolated positive eigenvalues of finite multiplicity and may also 
contain zero as an eigenvalue of possibly infinite multiplicity. Let s > 0. Given a 
function g : (7{E) —>■ C we define 

|^7(E)U,= ( J2 \9{K.[EW>^n[EYY\l<p< 00, 

A„[E]^0 

\g{E)\^^s= sup \g{Xn[E])\, 

where, as usual, each positive eigenvalue is repeated according to its multiplicity. 

The next lemma involves the Schatten norms ||-||c'-;l^''"<oo. For a compact 
operator £^ on a Hilbert space they are defined by H-Ellc = {J2n f^^iEYY^'^ , if 
r < 00, and H-EHcoo = ll-^'ll) where Hn{E) are the singular values of E, i.e., the 
non-zero eigenvalues of {E*EY^'^; recall that the Schatten space C, defined as the 
space of those compact operators for which the Schatten norm || ■ is finite, is a 
Banach space; see Reed and Simon [2Tj or Simon ^23J| for details. 

Let F := TT*. Recall that a{F) \ {0} = a{H) \ {0}, see [E]. In the next 
lemma, g{H) and g{F) are operators defined in the standard way by functional 
calculus. The following lemma is a variant of Lemma 8 of [2] . 
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Lemma 7 Let qo > 2, ■y > 0, p > qo/{qQ — 2), 2 < r < oo and s = 2qQj /[p{qo — 2)]. 
Then the following statements hold: 

(i) If the eigenfunctions of H satisfy (PI) then for any measurable function 
R : Q C and any function g : o"(if) ^ C we have 

\\Rg{H)\\cr < \\R\\L.r^n){\nr^\gm + C^^\g{H)\r,s). (4.5) 

(a) If the eigenfunctions of H satisfy (P2) then for any measurable matrix-valued 
function R on Q and any function g : cr{F) C such that if ^ cr(-^) then 
g{0) = 0, we have 



\\Rg{F)\\cr < C^^^Ml^^^^ \\R\\L.r^n)\9{F)\r,s. (4.6) 

Proof. We only prove statement (ii) since the proof of (i) is simpler. It is enough to 
consider the case where R is bounded and g has finite support: the general case will 
then follow by approximating i? in || ■ Hipr^fj) by a sequence n G N, of bounded 
matrix-valued functions and g in \ ■ |r,s by a sequence n G N, of functions with 
finite support, and observing that by (14.61) the sequence i?„(?„(F), n G N, is then 
a Cauchy sequence in C^. 

Since R is bounded and g has finite support Rg{F) is compact, hence inequality 
(14. 6 p is trivial for r = oo. Thus it is enough to prove (14. 6 p for r = 2 since the 
general case will then follow by interpolation (cf. [23]). It is easily seen that 

■= T'ilj^[H]/\\T,p^[H]\\ = Xn[HY^/^T^n[H], for all n G N such that X^[H] ^ 0, 
are orthonormal eigenfunctions of F, Fzn = \n[H]Zn, n G N, and span{2;„} = 
Ker(F)-'-. Hence, since 5^(0) = 0, 

oo 

\\R9m\h = Y.\\^3{F)z4l,^^) (4.7) 

n=l 

00 

= $^|^7(A„[i/])n|i?^.||i.(^) 

n=l 

00 

= Y.^rm-'\9{MH])?\\Ra'''yMH]\\l^in) 



n=l 



L2p/(p-i)(n) 



n=l 



00 

< C^||ai/lio.(n)||i?||i2.(^)^|^7(A4if])pA„[i/]^, 



n=l 



where for the last inequality we have used (14. 4p . This proves (14. 6 p for r = 2, thus 
completing the proof of the lemma. □ 

Recall that (5^(0, 0), 1 < p < 00, is defined in (II. 3p . 
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Theorem 8 (stability of resolvents) Let (A) be satisfied. Let ^ E C \ {cr{H) U 
{0}). Then the following statements hold: 

(i) There exists Ci > depending only on N, t, 6, a, c* and ^ such that ifSoo{(f), 0) 
< c^^, then ^ ^ o"(-ff) and 

II {w-'Hw -0''-{H- 0'' llc" < ci(5oo(0, 0). (4.8) 

(a) Let in addition (P) be satisfied by the operators H , H and T*ST for the same 
go, 7 and C. Let p > qo/{qo — 2) and r > max{2, a + ^^^^}- Then there 
exists C2 > depending only on N, r, 9, a, a*, r,p, go, C, 7, |f2| and ^ such that 
if ^pr{'P-i 4') < C2 ^, then ^ ^ cr(if) and 

\\{w'^Hw - 0"' -{H- i)-^\\cr < 02^.(0, 0). (4.9) 

Remark 9 Let s = [go/(go — 2)]max{2,a; + 27}. It follows by TheoremlE (H) 
(choosing p = go/(go — 2)j that if 5s{(l), 0) < c^"*^, then ^ ^ cr(-^) and 

Wiw-'Hw - - {H - < C2Ss{^,4>)- (4.10) 

Remark 10 As we shall see in Section 7, the best range for s in ^-lOj ) used in 
our applications is s > N; this corresponds to the case in which go = 00, 7 = N/A 
and a > N/2. See RemarksUi 

Proof of Theorem\^ In this proof we denote by c a positive constant depending 
only on A^, r, 6', a and c* the value of which may change along the proof; when 
dealing with statement (u) constant c may depend also on r, p, go,C, 7, \VL\. We 
divide the proof into two steps. 

Step 1. We assume first that ^ ^ cr(-^) U a{T*ST) and we set 

d,{i) = dist(^, a{H) U a{H) U a{T*ST)). 

In this first step we shall prove (14.81) and (14. 9 p without any smallness assumptions 
on 5oo(0, 0), (5pr(0, 0) respectively. 

We first prove (14. 8p . We shall use Lemma [3] and to do so we first estimate the 
terms Ai,A2, A^ in identity (14. ip . Clearly we have that 

<fi^ lel, V (4.11, 



|A„|i/]-^|-V d{i;,a(H)) 

Since the eigenvalues of the operator wT*STw coincide with the eigenvalues of H 
(see (Ol) ). it follows that 

00 ^ 



lel" ' ^V^ ' d{^,aiH))^ 
12 



^ +cfl + — ^L-l . (4.12) 



By (13. 6p and by observing that 

_ r _ >| N-l 

|detV0 - detV0| < N\N | V0 - V0| max | |V0| , | V0| | (4.13) 
it follows that 

|l-w|,|^^;-^^;~^| <c|V0- V0|. (4.14) 
Combining inequalities (I4.12p and (14.141) we obtain that 

\\Ai\\c«, WA-iWc < c (l + i| + ^^) ||V0 - V,^|U»|n) , 

We now estimate the term B in (14. ip . We recall that F = TT* and we set 
Fs = S^/'^TT* S^^"^ . Then, by polar decomposition, there exist partial isometrics 
Y,Ys : L'^{n,gdx) {L\Q,gdx))^ such that T = F^/'^Y and S^''^T = F^^Ys- 
We then have 

B = Y*Fg^\Fs - 0"'5^/'(^"' -I)iF- O'^F^/^Y . 

Hence, by Holder's inequality for the Schatten norms (see [211 P- 41]) it follows 
that 

||5||c" < ||f^/'(F5-0-'IIc-I|5'/'(5-^-/)IIloo(c)||(f-0"'^'/'IIc-. (4.16) 

Since ct(F) \ {0} = cr{H) \ {0}, we may argue as before and obtain 



2a 



Now, one easily sees that 

< c(|V0- V0I + |Ao0- Ao0|). (4.18) 
Combining fHTTCl) . (HTTP and fHTTg]) we conclude that 

By Lemma [3] and estimates (I4.15P and (I4.19p , we deduce that 



silc- <<;| i + M'*,'^')- (4.19) 



(^.■-'Hu.-0-^-(H-0-%^ < ( 1 + ^ + + J|j ) SM (4.20) 
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We now prove (14. 9p . In order to estimate Ai, A2, A3 we use estimate (14. 5 p and 
we get 



l^il 



A2\\cr <C( 1+^ 



|V0- V0||lp. 



3 



< c 



1 + lep 



II V0 - V0||Lpr(Q) 



(4.21) 



We now estimate B. We shall assume without loss of generality that 5* ^— / > 0. 
Thus, in order to estimate the norm of 5, we shall estimate the C^^ norms of 
Fs^\Fs - 0~^S^/\S-^ - lyl^ and {S'^ - If'^{F - i)-^F^I\ By Lemma [7] it 
follows that 



\\{S-'-lf/\F-i)-'F''^%,. 

00 

<c||(5-^-/)^/2|iiV(n)E 



ra=l 



< c\\S-^-I 



iLP'-(n) 



2r 



2go7 
(90-2) " 



(4.22) 



The same estimate holds also for the operator F^''^(F5—,^) -"^ 5*^/^(5 — J)^/^. Thus 
by Holder inequality for the Schatten norms it follows that 



|S||c. < c 1 + 



\s 



-1 



(4.23) 



By Lemma [3] and combining estimates and (I4.18p . (I4.2ip and (I4.23p we deduce that 



w-'Hw-O' -{H-O^'Wcr < ci 1 + — 



1^1 d^iO d^iO' 



). (4.24) 



Step 2. We prove statement (i). First of all we prove that there exists c > 
such that if 



5. 



< 



c(l + |eP + ci(e,a(i/))2) 



then ^^a{H)U a{T*ST) and 



d(e,a(#)), t/(e,a(T*Sr))> 



rf(e,a(i/)) 



(4.26) 



We begin with T*ST. By recalling that B = {T*ST - - {T*T - 0"^ (see 
the proof of Lemma [3]), by estimate f l4.19p with ^ = — 1 and by inequality f l4.35p it 
follows that there exists Ci > such that for all n G N 



1 



K[T*ST] + l A„[iJ] + l 



(4.27) 
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Assume that n G N is such that Xn[T*ST] < \^\ + d{^,a{H)); by KTl} it 
follows that if 

Ci(l + lei + d{^, a{H)))6U<P, 4>) < ^ ^^^^^ 



2(|e| + rf(e,a(iJ))) + l 
then 



Xn[H] < 



l-Ci[l+|e|+t^(e,^^(^))]5oo(0,</)) 

< 2(1^1 +d{^,a{H))), (4.28) 

(the elementary inequality {A + t){l-ty^ < 2A ii < t < A{2A + iy^ was used). 
Thus by K27^ and fOSD it follows that if 



2C^[1 + lel + rf(e, a(/7))] [1 + 2(|^| + d(^, a(i7)))] 
then 

|e-A„[r*5T]| 

> \^~K[H]\-\K[H]-K[T*ST]\ 

> rf(e, a(if)) - + lel + ci(e, a(if))][l + 2(1^1 + d{^, a(if)))]5oo(0, 0) 

> ^ (4.29) 

for all n G N such that A„[T*^r] < |^| + d{C,(T{H)). Thus inequality (14:261) 
for d{^,a{T*ST)) follows by ( 14.29^ and by observing that if n G N is such that 
Xn[T*ST] > lel + d{^,a{H)) then 1^ - Xn[T*ST]\ > d{^,a{H)). Inequality Km 
for d{^, (j{H)) can be proved in the same way. Indeed it suffices to observe that by 
Step 1 there exists C2 > such that 

<C2(5oo(0,0); (4.30) 



Xn[H] + l A„[i/] + l 

we then proceed exactly as above. 

By (14.201) and (I4.26P it follows that there exists c> such that if 

'-^^'^^-c(i + iep + .(e,a(i/))T ^'-''^ 

then i i a{H) U (7{T*ST) and 

(4.32) 

This completes the proof of statement (i). 

The argument above works word by word also for the proof of statement (ii), 
provided that (5oo(0, 0) is replaced by Spr{4>, 0). □ 
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Remark 11 The proof of Theorem\^ gives some information about the dependence 
of the constants ci,C2 on ^ which will he useful in the sequel. For instance, in the 
case of statement (i), in fact it was proved that there exists c depending only on 
N, T, 6, a, c* such that if ( [^^ . 31\ ) holds then ^4-32 ) holds. Exactly the same holds for 
statement (ii) where a depends also on r,p,qo,C,'-y,\Q\. Moreover, for such (f),(j), 
if ^ ct{H) then ^ and the summand 1/|^| + l/d{^,a{Q)) can be removed 

from the right-hand side of ^.32] ); furthermore, in this case statements (i) and (ii) 
also hold for = 0. This can be easily seen by looking closely at the proofs of (JjJ^ 
and fil21\ ). 



Remark 12 By the proof of Theorem\E it also follows that for a fixed ^ G C\ [0, oo[ 
no smallness conditions on Soo{4>, (p), Spr{(j), (p) are required for the validity of fi4-^ ) 
Ili4.9^ respectively. 



Theorem 13 (stability of eigenvalues) Let (A) be satisfied. Then the following 
statements hold: 



(i) There exists Ci > depending only on N, r, 6, a and c* such that ifSoo{4>, 0) ^ 



^, then 



E 

I n=l 



Xn[L] + 1 Xn[L] + 1 



l/o 



< Ci5a 



(4.33) 



(ii) Let in addition (P) be satisfied by the operators L, L and L for the same 
go > 2, 7 > and C > 0. Let p > go/(go - 2) and r > max{2, a + ^^^}- 
Then there exists C2 > depending only on N,T,6,a,c*,r,p,qo,C,'-f and \Q\ 



such that if 6. 



pr \ 



< Co 



then 



E 



vn=l 



l/r 



\n[L] 



K\L] + 1 



< C2(5 



pr \ 



(4.34) 



Proof. The theorem follows by Theorem [8] and by applying the inequality 

l/r 



E 



with El = vj-^Hw, E2 = H (see [231 P- 20]). 



< \\iEi + l)-'-{E2 + l)-%r, (4.35) 



□ 



Remark 14 We note that in the case of Dirichlet boundary conditions, i.e., V 



Wq' (fi), inequality li4-33\ ) directly follows by /5, Lemma 6.1] the proof of which 
based on the Min-Max Principle. 



IS 
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5 Stability estimates for eigenfunctions 



Definition 15 Let E be a non-negative self-adoint operator with compact resolvent 
on a Hilhert space H. Given a finite subset G of N we denote by Pg{E) the or- 
thogonal projector from Ti onto the linear space generated by all the eigenfunctions 
corresponding to the eigenvalues \k[E] with k E G. 

Observe that the dimension of the range of Pg{E) coincides with the number 
of elements of G if and only if no eigenvalue with index in G coincides with an 
eigenvalue with index in N \ G; this will always be the case in what follows. 

In the following statements it is understood that whenever n = 1 the term A^-i 
has to be dropped. 

Theorem 16 Let (A) be satisfied. Let X be a non-zero eigenvalue of H of mul- 
tiplicity m, let n E N be such that X = Xn[H] = ... = Xn+m~i[H], and let 
G = {n, n + 1, . . . , n + m — 1}. Then the following statements hold: 

(i) There exists ci > depending only on N, r, 9, a, c*, Xn-i[H], X, Xn+m[H] such 
that i/5oo(0, 0) < cj~^, then dimRanPc {w~^Hw) = m and 

\\Pg{H)-Pg{w-'Hw)\\ < ci5oo(0,0). (5.1) 

(a) Let in addition (P) be satisfied by the operators H , H and T*ST for the same 
qo, 7 and C . Let s = [qo/{qo — 2)] max{2, a + 27}. Then there exists C2 > 
depending only on N,T,9,a,c*,qo,G,'y,\Q\ A„_i[if],A, Xn+m[H] such that if 
Ss{^,<P) < then dimRa,nPG{w~^Hw) = m and 

\\Pg{H)-Pg{w~^Hw)\\ < C26s{<P,4>). (5.2) 



Proof We set p = idist(A, {a{H) U {0}) \ {A}) and A* = A if A is the first non-zero 
eigenvalue of H, and A* = A„_i[if] otherwise. 
By Theorem [13] (i) it follows that 

\Xk[H] - Xk[H]\ < c{Xk[H] + l){Xk[H] + l)5oo(0, 0) . (5.3) 

This implies that there exists c > such that if 5oo(0,0) < c~^^k[H]/{Xk[H] + 1)^, 
then Xk[H] < 2Xk[H]. This together with (15. 3p implies the existence of c > such 
that if 5oo(0,0) < c-'mm{p,Xk[H]}/{Xk[H] + if then \Xk[H] - Xk[H]\ < p/2. 
Applying this inequality for /c = n — 1, . . . , n + m, we deduce that if 



c(A„+,„|J/l + 1)= 

then 



|A4^] - X\ < p/2 , y k E G 

\Xk[H]- X\>3p/2 , y k eN\G . (5.4) 
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Hence dimRanPG'(w ^Hw) = m and by the well-known Riesz formula we have 
that 

PG[H] = -^J^{H-0-'dC, (5.5) 
Pg[w-'Hw] = j^{w-'Hw - O-'d^ , (5.6) 

where r(^) = A + pe^^ 0<e <27i. Hence 

\\Pg[H]- Pg[w~^Hw]\\ < p snp \\{w-^ Hw - {H - ^)-^\\. (5.7) 

Let Ci be as in Theorem [8] (i). By Theorem [8] (i) and Remark [TT] and by observing 
that A - p < 1^1 < A + p and 1/|^| < 1/p for all ^ G T, it follows that if 

*»'*-^)S,(l + A5l|//|+p=) 

then ^ 

\\{w-'Hw - ^)-' -{H- 0"! < ci (^1 + i + 5oo(0, 0). (5.8) 

The proof of statement (i) then follows by combining (15. 7p and (15.81) . The proof 
of statement {ii) is similar. □ 



Remark 17 The proof of Theorem [T^ gives some information about the depen- 
dence of the constants ci,C2 on Xn-i[H], X, Xn+m[H] which will be useful in the 
sequel. For instance in the case of statement (i) in fact we have proved that there 
exists c > depending only on N, r, 6, a, c* such that if 

A mm{p, A*} 

'-^^'^^-c(l + p^ + A„,.m 

then ^ 

\\Pg{H) - Pg{w-'Hw) \\<c(^1 + p + j^ 6M, 0)- (5.9) 

Exactly the same is true for statement (ii) where c depends also on qQ,C,'-f, 

We are going to apply the stability estimates of Theorem [16] to obtain stability 
estimates for eigenf unctions. For this we shall need the following lemma. 

Lemma 18 (selection lemma) Let U and V be finite dimensional subspaces of 
a Hilbert space Ti, dim U = dim V = m, and letui,..., Um be an orthonormal basis 
of U . Then there exists an orthonormal basis Vi, . . . ,Vm of V such that 

\\uk - VkW < 5''\\Pu - PvW , k = l,...,m, (5.10) 
where Pjj, Py are the orthogonal projectors onto U, V respectively. 
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Proof. Step 1. Clearly \\Pu - Pv\\ < 2. If 1 < \\Pu - Pv\\ < 2 then estimate 
(15.101) obviously holds for any choice of an orthonormal basis f i, . . . , fm of so we 
assume that \\Pu — Pv\\ < 1- Let u & U, \\u\\ = 1. Then 

\\Pvu\\ = \\u + {Pv - Pu)u\\ > 1 - \\Pu - PvW > . (5.11) 

Letting z = Pvu/\\Pvu\\ we have ||2;|| = 1 and 



hence 

II 7/ II ^ 

\\Pu -Pv\\'> \\{Pu - PvM' = \\uf - WPvuf > " : 



and therefore 

\\u- z\\ < V2\\Pu - Pvl (5.12) 
Step 2. Assume that ||Pc/ — Py|| < 1/6, and 

We shall prove that 

K^fe,^/)| <3||Pt;- Py II, kj = l,...,m, ky^l. (5.13) 
Indeed, we have for k ^ I 

\{PvUk,PvUl)\ = \{PvUk-Uk,PvUl) + {Uk,Ui) + {Uk,PvUl-Ui)\ 
= I {{Pv - PuW, PvUl) + (Mfc, {Py - Pu)ui) I 

< 2\\Pu-Pvl 

and the claim is proved by recalling (15. lip . 

Step 3. One can easily see that since ||Pc/ — Py|| < 1, the vectors zi,...,Zm 
are linearly independent. Thus we can apply the Gram-Schmidt orthogonalization 
procedure, i.e., define 

Zk-T!lZl {zk,vi)vi 
vi=zi, Vk = - — , /c = 2,...,m. 

Note that for /c = 2, . . . , m, 

- YaZI {zk, vi)vi\\ ) \\zk- YaZI {zk, vi)vi\\ 
and 

k-l k-1 

Zk-^{Zk,Vi)vi >l-^\{Zk,Vi)\. 



1 > 



1=1 1=1 
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Hence if 

Y,\{zk,vi)\<l (5.14) 
1=1 

then 

H-'4< '^^J};;'\ (5.15) 

Also for s = k, ... ,171, 



{Z s,Zk) \ + E?=l^ \{Zk,V l) I 
1 - Ef=/ 



S'ie^ We shall prove that for all = 2, . . . , m 

\\vk-z4<3-5'-'\\Pu-Pvl (5.17) 
\{zs,Vk)\<3-5''-'\\Pu-Pvl s = k + l,...,m, (5.18) 

provided that 

\\Pu-Pv\\ < (5.19) 

We prove this by induction. If A; = 2 then by fl5.13p and fl5.19p \{z2,Vi)\ = 
\{z2,z,)\<3\\Pu-Pv\\ < I hence by (El5]), 

and by (I5.16P and (I5.13P for s = 3, . . . , m also 

i — 3||/1/ — iVIl 

Let 2 < k < m—1. Assume that inequalities fl5.17p and (15.180 under assumption 
(15.191) are satisfied for all 2 < j < k. By assuming the validity of (I5.19P for k + 1, 
by (15.151) we obtain 



\\Vk+l - Zk+l\\ < 



l-3(EL5^-i)||Pt/-Pv 



Since ^fL^ 5'"^ < 574, by fICTD with k replaced by + 1, 3 J^Li 5'"i^i/ - 
PvW < 1/2 hence \\vk+i - Zk+i\\ < 3 • 5'=||Pc7 - Py||. Similarly, by fl5T6|) and (15^131) 
for all s = k + 2, . . .171 



3\\P u-Pv\\+3{EU^'-')\\Pu - Pv 
l-3(E-=i5^"i)||Pt/-Pv 



< ^ . t^'T „^ " „ < 3 ■ 5'\\Pu - Pv 
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Step 5. To complete the proof, we note that by (15.121) we have \\Ui — vi\\ < 
V2\\Pu-Pv\\- For A; > 2 we have that if (I5A9D holds then 

\\uk - Vk\\ < \\uk - ZkW + \\zk - Vk\\ < {V2 + 3 ■ 5^"^)||Pc/ - Pv\\ , 

while if (15.191) does not hold then \\Pu — Pv\\ > 10/(3 ■ 5^) and therefore 

\\uk - VkW < 2 < 3 • 5''-'^\\Pu - PvW- 

This completes the proof of the lemma. □ 

Lemma 19 Let (A) be satisfied. Let X be a non-zero eigenvalue of H of multiplicity 
m and let n & N be such that A = A„[-?/] = . . . = Xn+m-i[H]. Then the following 
statements hold: 

(i) There exists Ci > depending only on N, r, 9, a, c*, A„_i[if], A, Xn+m[H] such 
that the following is true: if S 00 {(f), <j)) < cj~^ and ipn[H], . . . ,ipn+m-i[H] are 
orthonormal eig en functions of H in L'^{fl, gdx) , then there exist orthonormal 
eigenfunctions ipniH], ■ ■ ■ ji'n+m-ilH] of H in L'^{Q,gdx) such that 

\\MH]-MH]\\LHn) < ci5oo(0,0), (5.20) 
for all k = n, . . . , n + m — 1 . 

(a) Let in addition (P) be satisfied by the operators H , H and T*ST for the 
same qo, 7 and C . Let s = [qo/iQo — 2)] max{2, a + 27}. Then there exists 
C2 > depending only on N,T,9,a,c*,qo,C,'y,\Q\, Xn-i[H], X, Xn+m[H] such 
that the following is true: if 5s{4>, 4') < c^"*^ and il)n[H], ... , il)n+m-i[H] are 
orthonormal eigenfunctions of H in L'^{Q, gdx) , then there exist orthonormal 
eigenfunctions ipniH], . . . ,ipn+m-i[H] of H in L'^{Q,gdx) such that 

\\MH]-MH]\\LHn) < C25s(0,0), (5.21) 
for all k = n, . . . , n + m — 1 . 

Proof. We shall only prove statement (ii) since the proof of statement (i) is 
similar. We first note that (pk '■= w'^ipklH], k = n, . . . ,n + m — 1, are or- 
thonormal eigenfunctions in L'^{Q, g dx) of w~^Hw corresponding to the eigen- 
values Xn[H], . . . , Xn+m-i[H]. By Theorem [16] and Lemma [18] there exists c > 
such that if Ss{4>, 0) < then there exist eigenfunctions ipn[H], . . . , ipn+m-ilH] of 
H corresponding to the eigenvalue A such that 

\\MH]-Vk\\LHn)<c6,i(l),4>). (5.22) 
In order to complete the proof it is enough to observe that 

\\(pk-ipk[H]\\L2^^^ < ||1 - u^"lL-(n)||^fc[^^]||LW(=-2)(n) < c||V0 - V0||L=(n) . 

□ 

In the following theorem we estimate the deviation of the eigenfunctions ipk[L] 
of L from the eigenfunctions ipk[L] of L. We adopt the convention that ipk[L] and 
ipk[L] are extended by zero outside 0(fi) and 0(fi) respectively. 
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Theorem 20 (stability of eigenfunctions) Let (A) be satisfied. Let X be a non- 
zero eigenvalue of L of multiplicity m and /ei ri G N be such that A = Xn[L] = 
. . . = Xn+m-i[L]. Then the following statements hold: 

(i) There exists Ci > depending only on A^, r, 6^, a, c*, A„_i[L], A, A„+m[L] such 
that the following is true: if 5oo[(t>)4') < c]"^ and ipn[L]i ■ ■ ■ iipn+m-i[L] are 
orthonormal eigenfunctions of L in L^(0(f2)), then there exist orthonormal 
eigenfunctions ipn[L\., ■ ■ ■ ,ipn+m-i[L\ of L in L'^{(j){Q)) such that 

WML] - ML]\\mHn)u,p{Q)) ^ c(5oo(0,0)+ 

+ \\ML] o (p - ML] o ^mn) + WML] o - o 0||i.(f,)),(5.23) 

for all k = n, . . . , n + m — 1 . 

(a) Let in addition (P) be satisfied by the operators L, L and L for the same 
qo, 7 and C . Let s = [qo/{qo — 2)] max{2, a + 27}. Then there exists C2 > 
depending only on N, r, 9, a, c*, go, C, 7, |^|, A„_i[L], A, Xn+m[L], such that the 
following is true: ifSsifp, 0) < c]~^ andipn[L], . . . , ipn+m-i[L] are orthonormal 
eigenfunctions of L in L'^{(j){Q)), then there exist orthonormal eigenfunctions 
ijjn[L], . . . y'lpn+m-ilL] of L in such that 

WML] - ML]WLH4,{n)u4>m < c(5,(0,0) + 

+ WML] o - ML] o 4>WL\n) + WML] o<j)-ML]o 0||L2(f,)),(5.24) 

for all k = n, . . . ,n + m — 1. 

Remark 21 We note that if in addition the semigroup e^^^ is ultracontractive 
then the eigenfunctions are bounded hence 

WML] o <j) - ML] o Hmn) + WML] o <P - ML] o 4>Wmn) < c{X)\V\'/', 

where V = {x & Q : 7^ 

Proof of Theorem We set 

ipkiH] = ipklL] o 0, 

for all k = n, . . . ,n + m — 1, so that ipn[H], . . . , ipn+m-i[H] are orthonormal 
eigenfunctions in L'^{Q,gdx) of the operator H corresponding to the eigenvalues 
Xn[H], . . . , A„+m_i[if]. By Lemma [19] (i) it follows that there exists Ci > such 
that if (5oo(0, 0) < then there exist orthonormal eigenfunctions ipn[H], . . . , 
ipn+m-ilH] in L'^{fl,gdx) of H corresponding to the eigenvalue A such that in- 
equality fl5.20p is satisfied. We now set 

ML]=MH]o 
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for all k = n, . . . , n + m — 1 , so that ipn[L], . . . , ipn+m-i[L] are orthonormal eigen- 
functions in L^(0(f2)) of L corresponding to the eigenvalue A. Then by changing 
variables in integrals we obtain 

\\ML]-ML]\\m4>{n)) < \\ML]°4>-ML]°4>\\LHn) 

< c(^\\iJk[L] o - i)k[L] o 0||L2(f^) + WiJklL] ocj)-iJk[L]o (Ph^m^ 

= c(^\\^k[L]o4)-^k[L]o(j)\\L2(n) + \\MH]-MH]\\LHn))- 
In the same way 

WML] - ML]\\LH<t,m < c{\\ML] o 4> - ML] o (Phnn) + WMH] - MH]\\mn))- 
Hence estimates flCTjl . (15:21) follow by flOOj) . fICTD respectively. □ 

6 On regularity of eigenfunctions 

In this section we obtain sufficient conditions for the validity of conditions (PI) 
and (P2). We begin by recalling the following known result based on the notion 
of ultracontractivity which guarantees the validity of property (PI) under rather 
general assumptions, namely under the assumption that a Sobolev-type Embedding 
Theorem holds for the space V. 

Lemma 22 Let Q be a domain in of finite measure and V a closed subspace 
of W^''^{^1) containing Wq''^{^1). Assume that there exist p > 2, D > such that 

\\u\\Lp(n) < L)\\u\\wi,2(n), (6.1) 
for all u eV . Then the following statements hold: 
(i) Condition /13.B) is satisfied for any a > 

(a) The eigenfunctions of the operators H, H and T*ST satisfy (PI) with go = 
oo, 7 = 2{p-2) ' ''^here C depends only on p, D, r, 9, c*. 

Proof. For the proof of statement [i) we refer to [31 Thm. 7] where the case 
V = W^''^{Q) is considered. The proof works word by word also in the general 
case. The proof of statement (ii) is as in [3, Thm. 7] where it is proved that for 
the Neumann Laplacian property (PI) is satisfied if (16. ip holds: this proof can be 
easily adapted to the operators H, H and T*ST. □ 

We now give conditions for the validity of property (P2). We consider first the 
case when an a priori estimate holds for the operators L, L, which is typically the 
case of sufficiently smooth open sets and coefficients. Then we consider a more 
general situation based on an approach which goes back to Meyers jTS] . 
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The regular case 



Recall that an open set in satisfies the interior cone condition with the 
parameters R > and /i > if for all x E Q there exists a cone C Q with the 
point vertex congruent to the cone 

K{R,h) = |x G : < (E^')'^' < ^ < 

In this paper the cone condition is used in order to guarantee the validity of 
the standard Sobolev embedding. 

The next theorem is a simplified version of P, Theorem 5.1]. 

Theorem 23 Let R > 0, h > 0. Let U be an open set in R^ satisfying the interior 
cone condition with the parameters R and h, and let E he an operator in L'^{U) 
satisfying the following a priori estimate: 

there exists B > such that if2<p<N + 2 and if u E Dom(£') and 
Eu e LP{U), then u G W'^^'p{V) and 

||M||vy2.P(t/) < B {^Eu\lv{jj) + ||M||i2([7)) . (6.2) 

Assume that Eip = Xip for some A G C. Then there exists c > 0, depending 
only on R, h, N and B, such that for /i = 0, 1, 

< c(l + |A|)^+^||V^|U2(^). (6.3) 

Theorem 24 Let (A) he satisfied and let 0(i7) and 0(^^) he open sets satisfying 
the interior cone condition with the same parameters R, h. If the operators L, L 
satisfy the a priori estimate lid. ^) with the same B, then the operators H , H , T*ST 
satisfy property (P) with go = oo, 7 = N/A and C depending only on T,R,h,c*,6 
and B. 

Proof. Recall that H, H and T*ST are the operators obtained by pulling-back 
to Q the operators L, L and L respectively. Clearly L also satisfies the a priori 
estimate 06.21) . Thus, by Theorem [23] the eigenfunctions of the operators L,L,L 
satisfy condition (16.31) hence, by pulling such eigenfunctions back to Q it follows 
that the eigenfunctions of if, H, T*ST satisfy (PI) and (P2) with go = 00, 7 = N/4 
and C as in the statement. □ 

The general case 

Here we shall assume that V = clvKi.2(n)Vo where Vo is a space of functions 
defined on fl such that C^{^) C Vo C Vr^'°°(fi). Moreover, for aU 1 < g < 00 we 
set 

Vq = cVi,q(Q)Vo. 
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Let — Aq : Vq (Vq')' be the operator defined by 

{-AqU,ip) = I Vu- Vipdx, 



Jn 

for all u E Vq, E Vq'. 

The following theorem is a variant of a result of Groger [13]; see also [2]. 

Theorem 25 Let (A) be satisfied. Assume that there exists qi > 2 such that the 
operator I — Aq : Vq ^ (H')' ^'^^ ^ hounded inverse for all 2 < q < qi. Then 
there exist qo > 2 and c > 0, depending only on Vq, t and 9 such that if u is 
an eigenf unction of one of the operators H , H , T*ST and A is the corresponding 
eigenvalue then 

\\\/u\\q < 0(1 + X)\\u\\q, (6.4) 

for all 2 < q < qo. 

Moreover, if Q is such that the interior cone condition holds then there exists 
c > depending only on V, r and 9 such that 

IIVmL <c(1 + A)||m||jv^, (6.5) 

N+q 

for all 2 < q < go- 

Proof. We prove the statement for the operator T*ST, the other cases being 
similar. We divide the proof into three steps. 
Step 1. We define 



Q{u,il>) = I uipgdx+ I aVu-Vipgdx, 
Jn Jn 



Qo{u,ip) = / uTpdx + I Wu-Wipdx, 
Jn Jn 

for all u e Vq, ip E Vq>. Sinctll 

\Qo{u,ilj)-f3Q{u,-^)\ < max{||l-/3^||i^(f^), ||/ -/3a^||L-(n)}||M|ki.9(n) 
there exist f3 > and < c < 1 depending only on A^, r and 9 such that 

\Qo{u,ip) - (3Q{u,ij)\ < c\\u\\wh<i^n)Mw^yin)^ (6-6) 

for all u E W^''^{n) and ^ E W^'^'i^). 

Step 2. Using the fact that ||(/ - As)"^ = 1, that q - Ag)~i|| is 

continuous and by observing that 2/(c + 1) > 1, it follows that there exists go > 2 
such that 

ll(/-A,)-^||<^, (6.7) 



2 Here we use ||/||^,.,(,,) = WfWl^^n) + II |V/| ||^.(^,) as the norm in W^'^^)- 
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for all 2 < g < go- By (16. 6p it then follows that for all 2 < g < go 

1 



inf 

ll^llwi 



sup Q{u,ip) > 



inf 



sup Qo{u,iIj)- 
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(6.^ 



Step 3. By (JHH) it follows that the operator / + {T*ST)g of Vg to V^, defined by 

(/+(T*Sr),M,7A) =Q(m,V^) (6.9) 
has a bounded inverse such that 

2(3 



[I+iT*ST)g)-'\\ = ( inf sup Q{u,ij)) < - 



(6.10) 



(6.11) 



for all il) '^Vqi. 

Now, if r2 satisfies the interior cone condition, then the standard Sobolev em- 
bedding holds. Thus, if g > 2 then q' <2 < N , hence Vq' is continuously embedded 



Then (16.41) follows by (16. 9p . (16.101) and by observing that 

Q{u,il)) = {1 + X) I uijjgdx^ 



Nq' 



into L^-i' (n). By (161TD we have 

ll^lkMo < (l + A)||(/ + (T*5T)j(-i)| 



sup 



wpg dx 



2/3 

< 1 (1 + A)||5'||loo(s^)||m|| Nq_ 

1 — c ^ ^ LT^(n) 



sup 



and the last supremum is finite due to the Sobolev embedding. 



Nq^ (^.12) 
□ 



Remark 26 IfQ satisfies the interior cone condition then inequality Ii6. 1\) is sat- 
isfied with p = 2N/{N — 2) if N > 3 and with any p > 2 if N = 2. Then by 
Lemmal2E it follows that condition liS. 2\) holds for any a > N/2 and the operators 
H, H, T*ST, L, L, L satisfy property (PI) with go = oo, 7 = N/A if N > 3 and any 
'-f > 1/2 if N = 2. In fact, if N = 2 property (PI) is also satisfied for 7 = 1/2; 
this follows by 07]. Thm. 2.4-4] ^^i^d /3, Lemma 10]. Thus by the second part of 
Theorem {2^ it follows that both properties (PI) and (P2) are satisfied for some 
go > 2 and 7 = A^(go - 2)/(4go) for any N >2. 

If Q is of class C^''^ (i.e., Vt is locally a subgraph of C^'^ functions) with < 
v < 1, then inequality / 1 6'. 1\) is satisfied with p = 2{N + u — 1) / {N — u — 1) , for any 
N > 2 (see also J^])- Thus Lemma implies that condition h3.^^) holds for any 
a > {N + u — 1)/ (2z/) and that the operators H, H, T*ST, L, L, L satisfy property 
(PI) with go = cx) and = {N + u - l)/(4z/). 
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7 Estimates via Lebesgue measure 



In this section we consider two general examples to which we apply the results of 
the previous sections in order to obtain stability estimates via the Lebesge measure. 

Let Aij G L°°(R^) be real- valued functions satisfying Aij = Aji for all i,j = 
1,. . . ,N and condition (12. 2p . Let Q he a bounded domain in of class 
and let F be an open subset of dQ with a Lipschitz boundary in dQ (see Defini- 
tion [27| below). We consider the eigenvalue problem with mixed Dirichlet-Neumann 
boundary conditions 



where u denotes the exterior unit normal to dQ. Observe that our analysis com- 
prehends the 'simpler' cases F = dQ (Dirichlet boundary conditions) or F = 
(Neumann boundary conditions), as well as all other cases where F is a connected 
component of dQ (the boundary of F in dQ is empty). See [13] for details. 

We denote by A„[f2, F] the sequence of the eigenvalues of problem (17. ip and by 
^/^^[fijF] a corresponding orthonormal system of eigenf unctions in L'^{Q). In this 
section we compare the eigenvalues and the eigenfunctions corresponding to open 
sets Q and Q and the associate portions of the boundaries F C dQ and F C dQ. To 
do so we shall think of f2 as a fixed reference domain and we shall apply the results 
of the previous sections to transformations and defined on Q, where (p = Id 
and is a suitably constructed bi-Lipschitz homeomorphism such that Q = 0(fi) 
and f = 0(F). 

Before doing so, we recall the weak formulation of problem (17.10 on Q. Given 
F C dQ we consider the space W^'^IQ) obtained by taking the closure of C^{Q) in 
W^''^{Q), where C^{Q) denotes the space of the functions in C°°{Q) which vanish 
in a neighborhood of F. Then the eigenvalues and eigenfunctions of problem (17.11) 
on Q are the eigenvalues and the eigenfunctions of the operator L associated with 
the sesquilinear form Ql defined on W := Wy'^{Q) as in (12. 3p . 

Definition 27 Let Q be a bounded open set in of class C^'^ and let F be an 
open subset of dQ. We say that F has a Lipschitz continuous boundary dT in dQ 
if for all X G 9F there exists an open neighborhood U of x in R^ and G ^{U) 
such that 

0(t/n(fiUF)) = {a; G R^: \x\ <1,Xn< 0}U{x G R^: \x\ < I, Xn < 0, Xi > 0} . 
7.1 Local perturbations 

In this section we consider open sets belonging to the following class. 

Definition 28 Let V be a bounded open cylinder, i.e., there exists a rotation R 
such that R{V) = Wx]a,b[, where W is a bounded convex open set in R^^-*^. Let 




on dn \ F, 



(7.1) 
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M, p > 0. We say that a bounded open se i C belongs to C^' {V,R,p) %f Q 
is of class (7™'^ (i.e., Q is locally a subgraph of C^'^ functions) and there exists a 
function g G C™'i(Vr) such that a+p <g<b, \g\m,i ■= J2o<\a\<m+i < 
M, and 

R{nnV) = {{x,xn) ■■ X eW , a < xn < g{x)}. (7.2) 



Let n,n e C°^\V,R,p) be such that n {VpY = Cln {VpY. We shall assume 
that the corresponding sets F C dQ, T C dCl, where Dirichlet boundary conditions 
are imposed, are such that 

rn\/'= = fn\/^ and p^(-i)^(rny) = p^(-i)^(f ny), (7.3) 

where -Pr(-i)vi/ denotes the orthogonal projector onto R^^^^W. Given F, condition 
(17. 3p uniquely determines T. 

Theorem 29 Let n G C°/(V, R, p) and let V be an open subset ofdVt with Lipschitz 
continuous boundary in dQ. Then there exists 2 < go ^ oo such that for any 
r > max{2, A^(go — 1)/Q'o} the following statements hold: 



(i) There exists ci > such that 



E 



1 



A„[f],f] + i A„[fi,r] + i 



l/r 



< Ci\VL A Vl\ -"0 , (7.4) 



for alip G Cl'j\V,R,p) such thatnn{VpY = nn{VpY, \n An\< where 
r C dCl is determined by condition ( [7.5] ]. 

(a) Let X[Q,r] be an eigenvalue of multiplicity m and let n & N be such that 
X[fl,T] = Xn[fi,T] = ■ ■ ■ = A„+m-i[^,r]- There exists C2 > such that the 

following is true: if n e c^^\v,R,p), n n {VpY = nn {VpY, \n An\< 

and r C dQ is determined by Iji7.3\ ) then, given orthonormal eigenf unctions 
ipn[^, r], • • • 5 i^n+m-il^, f]; thcrc cxist Corresponding orthonormal eigenfunc- 
tions 'j/'nffi, r], . . . , ipn+m-il^, T] such that 

_ _ _ qO-2 

\\M^,'^]-M^,'^]\\L^inun) < C2\n A fi| . 

Moreover, if in addition Aij G C°'^(R^), G C\^{y.,R., p) and V is a connected 
component of dQ then statements (i) and (ii) hold with Qq = oo. 

For the proof we need the following variant of jSj Lemma 4.1]. 

Lemma 30 Let W be a bounded convex open set in R^~^ and M > 0. Let < 

p < b — a and gi, g2 be Lipschitz continuous functions from W to H such that 

a + p < gi{x),g2{x) <b, (7.5) 
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for all X G W , and such that Lipgi, Lip5f2 < M . Let 6 = 2{b-a) '^^^ 3^ ~ 
min{5(i, 5(2} - 5|5fi - 5-21- Let 

Ok ■■= {{x,xn) ■■ X eW, a < xn < gk{x)} (7.6) 

for k = 1,2, 3. Let $ be the map from Oi into O2 defined as follows: 

if g2{x) < gi{x) then 

^(^'^^) - I f-^i-, ^ s , 'I S^'^^^l ^ S% 7n (7-7) 

\ {x,g2{x) + siiixN-gi{x))) if {x,xn) eOi\03, ^ ' 

while if g2{x) > gi{x) then 

Then ^ O3 C Ci n O2, 

\{xeOi: $(x)^x}| = |Oi\03|<2|Oi AO2I, (7.9) 

and $ zs a bi-Lipschitz homeomorphism of Oi onto Moreover $ G $t(^) 

where r depends only on N, M, 6. 

Proof. The proof is as in [Sl^ Lemma 4.1] where the case g2 < gi was considered: 
here we simply replace gi — g2 \gi — g2\- ^ 

Proof of Theorem We shall apply Theorems [13] and [20] with cf) = Id and 
given by 

X, X E VL\V, 



't'^'')-\R{-i)o^oR{x), xen^v. ^^-^^^ 

Here $ is defined as in Lemma [30] for gi = g and g2 = g, where g,g are the 
functions describing the boundary in V of Q, Q respectively, as in Definition [281 
Then clearly 4',4> E $r(^), where r depends only on N, V, M, p. Clearly 0(fi) = fl 
and 0(fi) = fi. Moreover, 0(r) = F, hence 

Moreover, condition (13.21) is satisfied for any a > N/2, see Remark [1] Hence 
assumption (A) is satisfied. Observe that by (17.91) and by the boundedness of the 
coefficients Aij, 

5p{(t),4)f <c I (|V0- V0|P + |Ao0-Ao0|P)dx < c|fi A (7.11) 

J {x<^Q.:(t)(x)^^{x)} 

By Theorem 3] the assumption of Theorem [25] is satisfied for the space Vq = 
C^{Vt) for some 2 < qi < 00. Thus by Remark [2^ the operators L, L and L 
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satisfy properties (PI) and (P2) for some 2 < go < oo and 7 = N{qQ — 2)/(4go)- 
Thus statement (i) follows by Theorem [13] (ii) with p = qo/{qo — 2). Moreover, 
Theorem [20] (ii) provides the existence of orthonormal eigenf unctions V^^ff^jr] 
satisfying estimate (15.241) with s = [go/ (q'o~"2)] max{2, A^(go— l)/go}. By Lemma [22] 
the functions V'fef^jr], V'fcl^jT] are bounded, hence by (17. 9p 

\\M^,^]°<I>-M^^^]°4>\\hin), \\M^,^]°<I>-M^,^]°^l^n) < A n\. (7.12) 

Thus statement (ii) follows by estimates (15.241) and (17.121) . 

Finally, if Aj e C^^\R^), Q,Q e C]}\V, R, p) and T is a connected component 
of dVL, by Troianiello Thm. 3.17 (ii)] the operators L and L satisfy the a priori 
estimate (16.21) on Vt and Vt respectively. Thus by Theorem [2l] the operators L, L 
and L satisfy properties (PI) and (P2) with go = 00 and 7 = A^/4, and the result 
follows as above. □ 



7.2 Global normal perturbations 

Let be a bounded domain with boundary. It is known that there exists 
t > such that for each x e {dVL)* := {x G : dist(a;, Sfi) < t] there exists 
a unique couple (x, s) e dVLx] — t,t[ such that x = x -\- suix). Clearly x is the 
(unique) nearest to x point of the boundary and s = dist{x, dfl). One can see 
that, by possibly reducing the value of t, the map x {x, s) is a bi-Lipschitz 
homeomorphism of {dQY onto dflx] — t,t[. Accordingly, we shall often use the 
coordinates {x, s) to represent the point x G {dQy. 

In this section we consider deformations of of the form 

Q = {Q\ {dQY) U {{x, s) G [dny : s < g{x)} (7.13) 

for appropriate functions g on dfl. 

Definition 31 Let Q and t be as above. Let < p < t and M > 0. We say that 
the domain Q belongs to the class Cj^'^{Q,t, p), m = or 1, ifVt is given by ( [7. j3[ j 
for some C™''^ function g on dQ which takes values in ] — t + p,t[ and satisfies 

\g\m,i < M. 

Given F C d^l and G C^'^{n,t, p), the set f C d^l where homogeneous 
Dirichlet boundary conditions are imposed, will be given by 

t = {{x,g{x)): XGF}. (7.14) 

Theorem 32 Let Q be an open set of class and t > be as above. Let F be 
an open subset of dQ with Lipschitz continuous boundary in dQ. Then there exists 
2 < go < 00 such that for any r > max{2, A^(go — l)/9o} the following statements 
hold: 
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(i) There exists ci > such that 

1 



E 



.n=l 



r^ V*^ 



~ 90-2 

<ci\nAn\^, (7.15) 



for all G C2}^(fi,t,p) snc/i i/iai, |r2 A i7| < c^^, where T C ^ is given by 



(a) Let A[f2,r] he an eigenvalue of multiplicity m and let n & N be such that 
X[fl,T] = Xn[fl,T] = ■ ■ ■ = A„+m_i[r2, r]. There exists C2 > such that the 
following is true: if VL E C^{Q,t, p), \fl A fl\ < €2^, and f C dQ is given 
by ^7.14^ then, given orthonormal eigenf unctions '0n[i^, T], . . . , ipn+m-i^, T], 
there exist orthonormal eigenf unctions F], . . . , ipn+m-i^i T] such that 

_ _ _ gO-2 

Moreover, if in addition Aij G C*^'"'^(R^), Q G Clf{Q,t,p) and T is a connected 
component of dQ then statements (i) and (ii) hold with Qq = 00. 

Proof. The proof is essentially a repetition of the proof of Theorem [29j the trans- 
formation $ is defined as in Lemma [30], with dQ replacing W and curvilinear 
coordinates (x, s) replacing the local euchdean coordinates {x,xn)- n 

8 Appendix 

In this section we briefly discuss how Theorem [8] can be used to obtain stability 
estimates for the solutions of the Poisson problem. 

Theorem 33 Let (A) be satisfied. Let the operators L, L, L satisfy (P) and fl 
satisfy the interior cone condition. Let f G L^(R^) and let v G W, "U be such 

that 

'{L + l)v = f, inm, .g^. 
iL + l)v = f, in0(fi). 

Let s = [go/ {qo — 2)] max{2, a + 27}. If N > 3, then there exists c > depending 
only on N,T,a,c*,qQ,C,'-y,fl such that 



h -v\\LHHn)uHn)) < ^ [m'^'' + Ss{<P, 0)) ||/||l2(r-) + 11/ o - / o j , 

where V = {x E Vl : 0(s) 7^ (/^ix)}. The same is true if N = 2 provided is 
replaced by \V\2~'' ^ e > 0. 

Proof. Observe that 

' {H + l){vo<p) = fo<p^ inn, 

{H + l){vo4)) = /o0, in ^ ' ^ 
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hence 

\\vo(l,-vo^L^n) < l|/o0-/o0||^.(f,) + ||(iJ+l)-i-(i/ + l)-i||||/o0|U.(f,). (8.3) 
By proceeding as in the proof of Theorem [20] one can easily see that 

11"" ~ '^llL2(0(Q)u0(f7)) 



< C \ \\V 0(j)~VO 0||l2(q) + \\vO(I)-VO 0||L2(f^) 

+ 11/ o - / o 011^2(0) + ||(^ + 1)-^ -{H+ o 0|U2(f,).) (8.4) 

By the Sobolev embedding it follows that if > 3 

\\vO(j)-VO (^||i2(f^), \\vO(t)-VO 0||i2(f^) 

< c|I?|l/^(||M|U2(n) + ||Vlz|U2(f,)) < c|P|l/^||/|U2(RiV). 

The same is true for N = 2 provided |X>|^/^ is replaced by l^l^^*^, e > 0. Moreover 
by Theorem [H] it follows that 

\\{H+ir' - (//+ i)-iii/o0ii^2(f,) < c5,(0,0)ii/iu2(K^). 

Thus, the statement follows by combining the estimates above. □ 

We now apply the previous theorem in order to estimate \\u — w||i^2(nuQ) where 

u, u are the solutions to the following mixed boundary valued problems and Vt is 
either a local perturbation of VL as in Section 17.11 or a global normal perturbation 
as in Section 17.21 

-Y.Z=^irXM^)t^)=f^ inn, 

M = 0, on r, (8.5) 

_d_ 

'i,j=l dxi 

u = 0, on r, (8.6) 

i,j=i ^tj ax, 



T.-,=.irXM^)§-) = f. inn, 



For any s > we set 

1/2 



Mf{s)= sup ( [ Ifl^dx] 



\A\<s 

The next theorem is a simple corollary of Theorem [33] and inequality (17. lip . 
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Theorem 34 Let Q, Q, T, T be either as in Theorem\2^ or as in TheoremlSR Then 
the following is true: there exists 2 < < oo such that for any r > max{2, N{qo — 
l)/go} there exists c > such that if \QAQ\ < then 

h - ^llL2(nuQ) <c[\nA n\ -0 ||/|U2(R^) + Mfic\n a . (s.r) 

Moreover, if in addition Aij G C^'^(R,'^), Q,Q E C^'^ and T is a connected compo- 



nent of dfl then estimate ( [g. 7| ) holds with qo = oo. 



Acknowledgments. This work was supported by the research project 'Prob- 
lemi di stabihta per operator! differenziali' of the University of Padova, Italy. The 
third author expresses his gratitude to the Department of Mathematics of the 
University of Athens for the kind hospitality during the preparation of this paper. 



References 

[1] G. Barbatis, Spectral stability under L^'-perturbation of the second-order 
coefficients, J. Differential Equations, 124, 1996, pp. 302-323. 

[2] G. Barbatis, Stability and regularity of higher order elliptic operators with 
measurable coefficients, J. London Math. Soc, 58, 1998, pp. 342-352. 

[3] V.l. Burenkov, E.B. Davies, Spectral stability of the Neumann Laplacian, 
J. Differential Equations, 186, pp. 485-508, 2002 

[4] V.I. Burenkov, P.D. Lamberti, Spectral stability of higher order uniformly 
elliptic operators, in Sobolev Spaces in Mathematics II. Applications in 
Analysis and Partial Differential Equations (to the centenary of Sergey 
Sobolev), edited by V. Maz'ya, International Mathematical Series, 9, 
Springer, New York, 2008, pp. 69-102. 

[5] V.I. Burenkov, P.D. Lamberti, Spectral stability of Dirichlet second order 
uniformly elliptic operators, J. Differential Equations, 244, pp. 1712-1740, 
2008. 

[6] V.I. Burenkov, P.D. Lamberti, Spectral stability of general non-negative 
self-adjoint operators with applications to Neumann- type operators., J. Dif- 
ferential Equations, 233, pp. 345-379, 2007. 

[7] V.I. Burenkov, P.D. Lamberti, M. Lanza de Cristoforis, Spectral stabil- 
ity of nonnegative selfadjoint operators. (Russian) Sovrem. Mat. Fundam. 
NapravL, 15, 2006, pp. 76-111. English translation in: Journal of Mathe- 
matical Sciences, 149, 2008, pp. 1417-1452. 

[8] V.I. Burenkov, M. Lanza de Cristoforis, Spectral Stability of the Robin 
Laplacian, Proceedings of the Steklov Institute of Mathematics, 260, 2008, 
pp. 68-89. 



33 



[9] E.B. Davies, Sharp boundary estimates for elliptic operators., Math. Proc. 
Cambridge Philos. Soc, 129, 2000, pp. 165-178. 

[10] E.B. Davies, Eigenvalue stability bounds via weighted Sobolev spaces., 
Math. Z., 214, 1993, pp. 357-371. 

[11] E.B. Davies, Heat kernels and spectral theory , Cambridge University Press, 
Cambridge, 1989. 

[12] P. Deift, Applications of a commutation formula, Duke Math. J., 45, 1978, 
pp. 267-309 

[13] K. Groger, A VT^'^-estimate for solutions to mixed boundary value prob- 
lems for second order elliptic differential equations. Math. Ann., 283, 1989, 
pp. 679-687. 

[14] J.K. Hale, Eigenvalues and perturbed domains. Ten mathematical essays 
on approximation in analysis and topology, pp. 95-123, Elsevier B. V., 
Amsterdam, 2005. 

[15] D. Henry, Perturbation of the boundary in boundary-value problems of 
partial differential equations, London Mathematical Society Lecture Note 
Series, 318, Cambridge University Press, Cambridge, 2005. 

[16] P.D. Lamberti, M. Lanza de Cristoforis, A global Lipschitz continuity re- 
sult for a domain dependent Dirichlet eigenvalue problem for the Laplace 
operator, Z. Anal. Anwendungen, 24, 2005, pp. 277-304. 

[17] P.D. Lamberti, M. Lanza de Cristoforis, A global Lipschitz continuity re- 
sult for a domain- dependent Neumann eigenvalue problem for the Laplace 
operator, J. Differential Equations, 216, 2005, pp. 109-133. 

[18] N.G. Meyers, An L^-estimate for the gradient of solutions of second or- 
der elliptic divergence equations, Ann. Scuola Norm Sup. Pisa, 17, 1963, 
pp. 189-206. 

[19] Y. Netrusov, Y. Safarov, Weyl asymptotic formula for the Laplacian on 
domains with rough boundaries. Comm. Math. Phys., 253, 2005, pp. 481- 
509. 

[20] M.M.H. Pang, Approximation of ground state eigenvalues and eigenfunc- 
tions of Dirichlet Laplacians. Bull. London Math. Soc., 29, 1997, pp. 720- 
730. 

[21] M. Reed, B. Simon, Methods of modern mathematical physics, II Fourier 
Analysis, Self-adjointness. Academic Press 1975. 

[22] G. Savare, G. Schimperna, Domain perturbations and estimates for the 
solutions of second order elliptic equations. J. Math. Pures Appl. (9), 81, 
2002, pp. 1071-1112. 



34 



[23] B. Simon, Trace ideals and their applications, Cambridge University Press, 
Cambridge- New York, 1979. 

[24] CM. Troianiello, Elliptic Differential Equations and Obstacle Problems, 
Plenum Press, New York, 1987. 



Gerassimos Barbatis 
Department of Mathematics 
University of Athens 
157 84 Athens 

Greece 

gbarbatis@math.uoa.gr 

Victor I. Burenkov 
Pier Domenico Lamberti 

Dipartimento di Matematica Pura ed Applicata 
Universitd degli Studi di Padova 

Via Trieste, 63 
35121 Padova 
Italy 

burenkov@math. unipd. it 
lamberti@math. unipd. it 



35 



